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We predict that stable mobile spatio-temporal solitons can exist in arrays of periodically curved 
optical waveguides. We find two-dimensional light bullets in one-dimensional arrays with harmonic 
waveguide bending and three-dimensional bullets in square lattices with helical waveguide bending 
using variational formalism. Stability of the light bullet solutions is confirmed by the direct numerical 
simulations which show that the light bullets can freely move across the curved arrays. This mobility 
property is a distinguishing characteristic compared to previously considered discrete light bullets 
which were trapped to a specific lattice site. These results suggest new possibilities for flexible 
spatio-temporal manipulation of optical pulses in photonic lattices. 
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I. INTRODUCTION 

Discrete solitons, the self-localized states which can 
exist in both one- and two-dimensional nonlinear pho- 
tonic lattices, have attracted a lot of attention over the 
recent years [l|. The interplay between the underly- 
ing periodicity and nonlinear self-action enables a whole 
lot of new opportunities for all-optical control, which 
are not possible in homogeneous media [1]. Of par- 
ticular interest are pulses which are localized both in 
space and in time as a result of balance between the 
self-focusing nonlinearity, diffraction, and anomalous dis- 
persion. Such spatio-temporal optical solitons are called 
light bullets [2j. Light bullets are inherently multi- 
dimensional, and consequently they are always unstable 
in Kerr media 0]. Discreteness in photonic lattices can 
provide an additional stabilization mechanism, and in- 
deed stable continuous-discrete solitons were found in 
nonlinear waveguide arrays [!, Hj]. However, these dis- 
crete light bullets are strongly trapped in a few waveguide 
channels, and hence their mobility is practically reduced 
to the longitudinal dimension, whereas they cannot move 
across the photonic lattice. 

Recently, it was suggested that arrays of periodically 
curved optical waveguides offer unique opportunities for 
the control of both the strength and frequency disper- 
sion of the waveguide coupling [B-Q • Propagation of light 
beams in various types of the nonlinear lattices with mod- 
ified diffraction has been analyzed 0-[3 • It was demon- 
strated that light beams can become localized in the peri- 
odically curved nonlinear waveguide arrays in the form of 
diffraction- managed solitons (lSMlTT l . These diffraction- 
managed spatial solitons are reminiscent of dispersion- 
managed temporal solitons (l8l - [20j . However, nonlinear 
propagation and localization of light pulses in diffraction 
managed photonic lattices have not been analyzed yet. 

In this work, we predict, for the first time to the best of 
our knowledge, that stable broad spatio-temporal solitons 
can exist in arrays of periodically curved optical waveg- 
uides, and such discrete light bullets are mobile and can 
move across the array. We find approximate light bullet 



solutions using the variational technique, where we as- 
sume a Gaussian shape of the spatial-temporal soliton as 
shown in Fig. H^b). The results of the variational anal- 
ysis indicate that by optimizing the waveguide bending 
amplitude and period, one can achieve stable and mobile 
two- and three-dimensional light bullets. We emphasize 
that we find stable light bullets which spatial width is 
much larger than the lattice period. Hence, the light 
bullets would not be trapped in a single waveguide, re- 
taining mobility in the transverse dimension. This con- 
clusion is confirmed by direct numerical simulations. We 
also demonstrate that a similar approach can be used 
to achieve stabilization of multidimensional matter-wave 
solitons in Bosc-Einstcin condensates. 



II. TWO-DIMENSIONAL LIGHT BULLETS 

A. Discrete-continuous model equation 

First, we consider propagation of light pulses in a one- 
dimensional array of coupled optical waveguides, where 
waveguide axes are periodically curved in the longitudi- 
nal propagation direction z, see the sketch in Fig. [TJa). 
We model the light bullet dynamics with the modified 
discrete-continuous nonlinear Schrodinger propagation 
equation which takes into the periodic waveguide bend- 
ing [12, [ll| and temporal dispersion [22hj25j : 
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Here n is the waveguide number, a n are the complex 
amplitudes of individual waveguide modes, z is the prop- 
agation distance, t is time, C is the coupling coefficient 
between the modes of straight neighboring waveguides, 
D is the second-order dispersion coefficient, 7 is the non- 
linear coefficient, and p(z) is the phase coefficient related 
to the waveguide bending [2l[. We consider all variables 
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in Eq. ([T]) to be dimensionless, applying the following 
normalization: 
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o« — A n / \/ J s , z — Zd/z s , t — td/t Sl 
C = C d z s , D = D d z s t~ 2 , 7 = j d I s z s , 



and finally obtain 
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where -A„ and variables with the subscript d denote the 
values in physical dimensions, and subscript s denotes the 
scaling coefficients. With these notations, the coupling 
phase coefficient is defined as 



We further consider the harmonic waveguide bending 
profile, 



xq(z) = xom cos(2nz/L), 



(7) 



p(z) = —— — x 
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where xo m is the bending amplitude (normalized to x s ) 
and L — L d /z s is the normalized bending period. To 
reduce the number of independent coefficients in Eq. ([6]) , 
we choose the scaling coefficients as follows: 



where xq(z) is the waveguide bending profile normalized 
to x s , prime stands for the derivative, d x — dd/x s is the 
normalized transverse separation between the neighbor- 
ing waveguides, Xd is the laser wavelength in vacuum, 
and no is the average medium refractive index. We note 
that the values of the scaling coefficients x s , z s , t s , and 
I s can be chosen arbitrarily. Accordingly, the results ob- 
tained for a specific set of dimensionless parameters can 
be related to a variety of physical configurations. 



Zs = *« = V\ D d\z s , x s = d d \f2C d z s , I s = 
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Then, the normalized coefficients are 

L = 2tt, d T — — — 
' 2C 2 



D = <td = sign(L>d) = ±1, 
7 = u 1 = sign(7d) = ±1, 
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B. 



Continuous approximation for mobile solitons and Ec l- © can be wri tten as 



As outlined in the introduction, our aim is to iden- 
tify the conditions for the existence of stable and mo- 
bile solutions in the form of light bullets, which remain 
localized in space and time as they propagate along the 
waveguide array. It was found that the solitons which are 
strongly localized in the spatial domain tend to become 
trapped at a particular lattice site [H, H3| • Only when 
the soliton width extends over five or more waveguides, 
the soliton becomes mobile 28]. However, such broad 
spatio-temporal solitons were found to be unstable [3, 0] . 

In order to study the key effects of periodic waveguide 
bending on the stability of broad spatio-temporal light 
bullets, we employ the continuous approximation. Specif- 
ically, we consider wavepackets that are slowly varying in 
space, and describe their profiles with a continuous func- 
tion, a n (t, z) = u(x = nd,t,z). We expand u in Taylor 
series in x and substitute in Eq. ([T]), 
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In self- focusing media we have ct 7 = 1, whereas ct 7 = — I 
in case of self-defocusing nonlinearity. In order to achieve 
soliton localization in temporal dimension we require 
that the temporal dispersion is either anomalous in self- 
focusing media , or dispersion is normal if the nonlin- 
earity is self-defocusing, i.e. ctdo^ > 0. 

We note that there is just one parameter, the cou- 
pling coefficient C, which defines the mapping of solu- 
tions based on the continuous Eq. (fTU| to the original 
discrete-continuous Eq. ([1]). Continuous solutions may 
serve as a good approximation provided that their spa- 
tial width (W) is much larger than the distance between 
the neighboring waveguides, 



+ j\ii\ 2 u = 0. 



2 dt 2 

Then, we make a transformation 

u{x, t, z) = u(x — Xx(z), z) exp 
where 



(4) 



w-»d x 



I 

2C 2 



2n 2 

r 2 r 2 



2L| 
LI 



(12) 



2iC I cosp(zi)dz\ 
o 



(5) 



x\{z) = 2Cd I sinp(zi)dzi, 
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where we introduce the notation L c — ir/C d to define 
the physical coupling length between two straight waveg- 
uides [20| . It follows that Eq. (fTOj) can be used as an ap- 
proximate model for Eq. (JT|) provided that the normalized 
coupling coefficient (C), or equivalently the ratio of the 
bending period to the coupling length, are sufficiently 
large. 
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FIG. 1. (a) Schematic of a periodically curved one- 
dimensional waveguide array, (b) Profile of a spatiotempo- 
ral light bullet with Gaussian ansatz. (c) Stability diagram 
in the two-dimensional case. The dashed lines indicate the 
value of the critical energy E CI and the boundary between 
positive and negative average diffraction regions (A), (d) The 
dependence of E vs. Atp (phase increment after one period) 
for B = 2.1. Solid line corresponds to stable solutions, dashed 
line to unstable solutions. The dash-dotted line indicates the 
diffraction map strength S. 



C. Variational approach 

We now look for approximate analytical solutions of 
the Eq. (fl"0|) using the variational method [30|. We con- 
sider the Gaussian Ansatz for the light bullet spatiotem- 
poral profile [see Fig.QJb)], 

1 " 
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u = A(z) exp |i0(z 

+ I[b(z)x 2 +I3(z)t 2 



W 2 (z) T 2 {z) 



(13) 



where A and <fi are the amplitude and phase of the soliton, 
T and W are its temporal and transverse spatial widths, 
and (3 and b are the temporal and spatial chirps. The 
Eq. (fTO)) can be derived from the Lagrangian 
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We substitute the ansatz Eq. (fT3")) into Eq. (fT4")l and in- 
tegrate over spatial coordinates to obtain the reduced 
Lagrangian 

(4/7r)£ cff = A 2 WT [4c// - b'W 2 - f3'T 2 

-\{z)b 2 W 2 - P 2 T 2 + A 2 /V2 



W~ 2 - DT~ 2 
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which variation with respect to <f> yields the energy con- 
servation dE/dz = with the constant E being propor- 



tional to the pulse energy 



E = A 2 WT 



\u\ dxdt. 
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The conservation of E is used to eliminate A from the 
remaining variational equations 



W = XbW, 
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where the equations for j3' and T' were used to obtain 
the single second order equation for T" . 

The systematic results obtained using the variational 
model Eq. (fT8)) are shown in Fig. [He). The dotted area 
in Fig. [ljc) corresponds to the values of parameters E 
and B for which stable light bullet solutions exist which 
oscillate with the period of the waveguide modulation. 
The solutions of the variational equations were obtained 
using Newton relaxation method. The two boundaries of 
the stability region can be found analytically, they are 
indicated in Fig. [ljc) by the two dashed lines . The 
first vertical line on the right corresponds to the value of 
waveguide bending B for which the average diffraction 
(A) is equal to zero. The bright soliton solutions ex- 
ist only for positive average diffraction since we assumed 
that the nonlinearity tr 7 is positive. When approaching 
this limiting value, the spatial extent of solutions grows 
towards infinity. The second limiting curve in Fig. [ljc) 
indicates the values corresponding to the critical energy 
E cr = 2yJ (A) of the Townes fundamental soliton [3ll - [33l | 
for the average diffraction (A). Stable solutions exist only 
for values of E > E CY . 

For a fixed value of the structural parameter B, there 
exists a family of soliton solutions. These solutions can 
be characterized by the pules energy E, the phase accu- 
mulated over one period A</> (this parameter has a similar 
meaning to the propagation constant for solitons in ho- 
mogeneous structures), and the diffraction map strength 
S. We define the values of S analogous to dispersion 
map str eng th in structures with varying temporal dis- 
persion [341 l35l| : 
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where WVwhm is the full width at half maximum of the 
wave packet at the point where it is the narrowest along 
the x direction. Typical dependencies of E and S on 
Acf> are shown in Fig. [Tfcl) . The value of E at A(p — > is 
equal to the critical energy of the Townes soliton E CI . We 
observe that the energy is a monotonously growing func- 
tion of A<ft, which indicates that the Vakhitov-Kolokolov 
criterion [36| for soliton stability is fulfilled. However this 
criterion is not sufficient to predict stability for period- 
ically modulated structures, and we find that the solu- 
tions become unstable for high values of the energy E 
and diffraction map S. 
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FIG. 2. Stable two-dimensional light bullet moving across the 
array. Parameters values are B — 2.1, C = 100, D = 7 = 1, 
and the initial Gaussian shape corresponds to the variational 
solutions with E = 1.1, W(0) = 0.68, T(0) = 1.55, 6(0) = 0, 
T'(0) = 0, multiplied by a spatial chirp function exp(ifcn) 
with k — 7r/20. (a) The spatial and temporal widths Wn 
and Tn, and the total energy En of the propagating pulse 
found by the numerical integration of the discrete propagation 
equation {TJ. (b)-(e) Snapshots of the spatiotemporal pulse 
shape at several propagation distances. 



D. Mobile bullets: full numerical modeling 



We note that the equation (pTJj) admits solutions mov- 
ing in both spatial and temporal dimensions. For exam- 
ple, any solution u(x, t, z) can be multiplied by a linear 
phase which gives a new solution u{x — X2) exp[ikx + 
i(f>2(z)], where 



dx 2 , , 
— = k\{z) 
az 



(21) 



Hence, pulses can move along ^-direction during the 
propagation, and this movement gives rise to an aver- 
aged soliton drift if the mean diffraction (A) is nonzero. 
However, the original discrete-continuous Eq. ([!} does 
not have such symmetry, and soliton motion can be sup- 
pressed due to the Peierls-Nabarro potential associated 
with the underlying discretness [26l428| |. Therefore, it is 
essential to confirm the soliton mobility through direct 
numerical simulations of Eq. ([I) . 

In Fig. [2] we present an example of a stable light bullet 
which moves across the array according to numerical sim- 
ulations of the full discrete-continuous propagation equa- 
tion ([T|). Here, the initial condition was taken from the 
periodic solution of the variational equations (|18[) multi- 
plied by a spatial chirp function exp(ifcn) with k = 7r/20, 



which corresponds to the initial beam tilt. As depicted in 
Fig. [H[a), the spatial and temporal widths of the soliton 
oscillate during the propagation in an irregular pattern, 
but dynamical stability is achieved over long propagation 
distances. The widths in the numerical simulations are 
calculated according to the formula 



(n}\: 



(22) 



where (n) = J2 n I n \ a n(t, z)\ 2 dt/ J2 n I \ a n(t, z)\ 2 dt. It 
gives approximately the full-width at half-maximum 
(FWHM) value for a Gaussian pulse. The numerically 
calculated energy is defined as 



En = d x 



E 



\a n (t,z)\ 2 dt. 
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Moving solitons in lattices are known to exhibit radia- 
tive losses 37]. In our simulations, we observe that the 
energy of the light bullet slowly decreases mainly due to 
radiation in the waveguide sections with the highest cur- 
vature. This radiation is more pronounced for tighter 
bending [smaller values of C corresponding to smaller 
bending period Ld acording to Eq. (fT2"j) ]. In Fig.[Hb)-(e) 
the spatiotemporal profiles of the light bullet are shown 
for several values of the propagation distance z. Here, the 
solution is shown in the reference frame moving along the 
periodic translation with zero mean shift resulting from 
the waveguide bending. The stable evolution presented in 
Fig. [2] corresponds to typical physical parameters for op- 
tical pulses Wd — 100 /iixi and Td — 50 fs in photonic 
lattice with the physical characteristics of dd = 5 /zm, 
Ld = 10cm, Cd — 6.3 mm -1 , and Dd = 175 fs 2 cm -1 . 

For comparison, in Fig. [3] we show a solution for the 
same values of the parameters as in Fig. [5] except for a 
lower pulse energy, corresponding to the unstable region 
in Fig. [He). I n this case, the pulse quickly spreads out 
in both space and time. 



III. THREE-DIMENSIONAL LIGHT BULLETS 

We now consider the possibility to stabilize three di- 
mensional light bullets, while preserving their mobil- 
ity. This is a more challenging problem, since two- 
dimensional light bullets in media with cubic nonlinearity 
exhibit critical collapse, and even weak perturbations can 
be sufficient to prevent the collapse effect [32j]. On the 
other hand, the three-dimensional light bullets can ex- 
hibit super-critical collapse, which can only be prevented 
through strong changes in the physical system 32 1. 

We consider the existence and dynamics of three- 
dimensional light bullets in the two-dimensional peri- 
odically curved square lattice of optical waveguides, as 
schematically shown in Fig. 0Ja). For the square lat- 
tice, the discrete-continuous propagation equation has 
the same form as Eq. ([!]) with the additional coupling 
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FIG. 3. An example of an unstable solution. Parameters are 
the same as in Fig. [2] except E — 0.6. 



along the y axis 

+ C (e-^a n+ltm + e*Wa„_ lim ) + 

+ C (e-^Vm+i + e^Vm-i) + (24) 
Dd 2 a„ 



2 dt 2 



+ l\a n ,m\ a n , m = . 



where y — yd/y s is the normalized coordinate, and we 
assume the same scaling of the transverse spatial dimen- 
sions, y s = x s . Then, the coupling phase coefficients are 
defined similar to Eq. ([3]), 

_ 2nn dx 2 , _ 2irn dx 2 

Px\z) — ^ x , p y \z) — - y , [zo) 



X d Z s 



X d z s 



where Xq(z) and yo(z) define the transverse bending pro- 
file. 

We further consider the helical waveguide bending pro- 
file [see Fig. 0Ja)], 

x (z) = x 0m cos(27rz/L), y (z) = y 0m sin(27rz/L), (26) 

where xo m = yo m is the bending amplitude (normalized 
to x s ) and L is the normalized bending period. 

Following the approach described above in Sec. Ill BL 
we consider wavepackets that are slowly varying in space 
and describe their profiles with a continuous function, 
0"n,m{t, z) — u(x — nd, y — md, t, z). Then, we apply the 
scaling coefficients according to Eq. (j9|), and derive the 
approximate continuous equation [cf. Eq. ()10|)]: 



du X x (z) d 2 u X v (z) d 2 u uu d 2 u 

(27) 



where \ x (z) and X y (z) the normalized diffraction coeffi- 
cients corresponding to the two transverse directions: 



X x (z) = cos(Bsinz), 
X y (z) = cos(B cos z), 
A-K 2 n dx 2 



(28) 
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Next, we perform the variational analysis using the 
three dimensional Gaussian Ansatz in the following form: 
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Following similar procedure as we used in the two- 
dimensional case [Sec. Ill C| . we arrive at the set of vari- 
ational equations 
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(30) 



E 



where E = A 2 WVT. 

We perform numerical analysis to identify periodic so- 
lutions of the variational equations (f30|) . Such solutions 
correspond to three-dimensional light bullets, where the 
light-bullet extension along all dimensions (x,y,t) is re- 
stored after each bending period. We present an example 
of such solution in Fig.[3Ib), which illustrates a stable and 
mobile three-dimensional light bullet. Indeed, the light 
bullet spatial and temporal localization widths are re- 
stored to the same values after each bending period. The 
solid line at the bottom shows the transverse position 
X2 of the pulse center as a function of the propagation 
distance z for a finite initial beam tilt. The light bullet 
wobbles in the transverse direction during the propaga- 
tion, and it is clear that on average the pulse moves across 
the array. 

We summarize the existence and stability properties of 
three-dimensional light bullets in Figs. E£c,d). We show 
the dependence between the energy E and the diffrac- 
tion map strength S on the phase increment A<f> for one 
bending period. In comparison with the two-dimensional 
case [Fig. (Hd)], this dependence is much more complex, 
forming a loop-like structure. Interestingly, stable so- 
lutions, depicted with a solid line, are present for both 
positive and negative curve slopes of the E(A<fi) depen- 
dence, which indicates that the Vakhitov-Kolokolov cri- 
terion [36| is not applicable here. 
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FIG. 4. (a) Sketch of helical waveguide array in three dimen- 
sions, (b) The dependence between the soliton energy E and 
the phase increment A(j> for one bending period. Solid line 
corresponds to the stable solutions, dashed line corresponds 
to the unstable solutions, (c) Stable three-dimensional light 
bullet moving across the helical array. The transverse spa- 
tial widths W (solid), V (dotted) and the temporal width T 
(dashed) of the light bullet are shown as functions of the prop- 
agation distance z. The solid line below shows the position of 
the moving pulse according to Eq. (|2ip for k x = 0.3. 



tive mass of the mater waves [38H40j . allowing for sta- 
bilization of dispersion-managed solitons in two dimen- 
sions [4l|. Stabilization of subdiffractive solitons with 
almost zero dispersion in lattices in two dimensions was 
also proposed [13( . Another approach consists of modula- 
tion of the nonlinear coupling using Feshbach resonances. 
In this context, stable two dimensional solitons were pre- 
dicted [42M44I . however stabilization in the full three- 
dimensional case requires an additional potential (45l - |47| , 
feedback control of the nonlinearity [48| or inclusion of 
inelastic losses [49]. Nonlocal dipolar interactions were 
also shown to provide stabilization of quasi two dimen- 
sional BEC solitons 0, HH. 

Our method can provide an alternative solution to the 
stabilization of matter-wave solitons. We estimate that 
in the case of a Bose- Einstein condensate of 7 Li atoms in 
the quasi two dimensional setting with the axial trapping 
frequency lu± — 2tt x 300 Hz, the lattice constant has the 
value d = A/2 = 0.4 /im, the width along the lattice axis 
W x — 6.5 jum, transverse width W y = 20.9 /im and the 
number of atoms N — 4000, where we assumed that the 
scattering length is tuned to a s = — 5ao using a Feschbach 
resonance [52j . The time period of the lattice modulation 
is equal to r = 45 ms and the total time of evolution in 
Fig. (Ha) is 0.7 s. 



CONCLUSIONS 



IV. MATTER- WAVE SOLITONS 

Finally, we note that a very related problem in the 
field of Bose-Einstein condensates is the stabilization of 
multidimensional matter- wave solitons, which dynamics 
can be described by the discrete GrossPitaevskii equation 
which is mathematically equivalent to Eq. ([1}. Several 
methods of stabilization have been considered to date. 
Discreteness can be introduced by applying deep opti- 
cal lattices, with the cost of reducing the soliton mo- 
bility to one dimension On the other hand, shal- 
low optical lattices can be used to manipulate the effec- 



We have demonstrated that arrays of curved waveg- 
uides can be used to stabilize two- and three-dimensional 
light bullets against the collapse in media with Kerr-type 
nonlinearity. Most remarkably, the bullets can freely 
move across the periodic structure, unlike the previously 
considered discrete light bullets which were trapped to a 
specific lattice site. These results suggest new possibil- 
ities for flexible spatio-temporal manipulation of optical 
pulses, and they are also relevant for the stabilization of 
multidimensional matter- wave solitons. 

This work was supported by the Australian Research 
Council. 
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